Unadulterated spectral function of low energy quasiparticles: Bi-2212, nodal direction 
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Fitting the momentum distribution photoemission spectra to the Voigt profile appears to be a 
robust procedure to purify the interaction effects from the experimental resolution. In application 
to Bi-2212 high-Tc cuprates, the procedure reveals the true scattering rate at low binding energies 
and temperatures, and, consequently, the true value of the elastic scattering. Reaching the minimal 
value ~ 16 meV, the elastic scattering does not reveal a systematic dependence on doping level, 
but is rather sensitive to impurity concentration, and can be explained by the forward scattering on 
out-of-plane impurities. The inelastic scattering is found to form well-defined quasiparticles with the 
scattering rate ~ lu 2 and ~ uj s , above and below T c , respectively. The observed uj — T asymmetry 
of the scattering rate is also discussed. 
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A century ago, in time of the renowned Kamerlingh 
Onnes's experiments, the search for actual behavior of 
electrical conductivity in pure metals at low temperature 
has led to discovery of superconductivity [l| . Today, the 
similar problem, the low temperature/energy behavior 
of normal electrical conductivity, is topical again, now 
for high-T c cuprates (HTSC). The quantity of special in- 
terest is a quasiparticle scattering rate 1/t @ or, more 
precisely, the single particle self-energy, E = E' + iE" [E], 
the real part of which can be associated with the mass 
renormalization and the imaginary part is proportional 
to the scattering rate or, in the simplest Drude model, 
to the normal state resistivity. The true lowest value of 
E" (taken at the Fermi level) is important to know in or- 
der to reconcile the parameters of quasiparticle spectrum 
with transport measurements, but its asymptotic behav- 
ior, i.e. £"(w,T) at low energy and temperature, is vital 
to judge whether the quasiparticle approach is applicable 
at all to describe the electronic properties of HTSC. 

The self-energy function, £(w,T), is closely related 
to (and can be derived from) the quasiparticle spec- 
trum, presented by the quasiparticle spectral function 
A(k, u), T), which, in turn, can be mapped accurately in 
the whole Brillouin zone (BZ) by modern angle resolved 
photoemission spectroscopy (ARPES) 4J. For the nodal 
direction, along which the superconducting d-wave gap 
function has a node (i.e. changes sign), all the interac- 
tions which form the nodal quasiparticles are encapsu- 
lated in E, both parts of which can be confidently deter- 
mined from ARPES spectra in the range about 0.3 eV 
below the Fermi level Nevertheless, the detailed be- 
havior of A(u>, T) and E(w, T) in the very vicinity to Ep 
remains puzzling. The experimental resolution, which 
can be safely neglected for higher binding energy, plays 
a crucial role here. Roughly, the photocurrcnt intensity 



can be well approximated by a convolution of the spec- 
tral function, multiplied by the Fermi-function f(u>), and 
overall experimental resolution 0: 



J(k, oj) oc A(k, u>)f(u) <g> i?(k, oj). 



(1) 



The latter consists of two essential components, the re- 
sponse function of the analyser, iiU(k, oj), and one which 
accounts for inhomogeneities of sample surface, i?s(k, lu): 
R = Ra <8> Rs ■ While the analyser response function is 
fixed and can be measured independently, the surface 
inhomogeneities (mechanical, chemical or in charge dis- 
tribution) result in a systematic error which is difficult 
to account for. This seemingly technical problem has 
rather strong fundamental impact, setting an unavoid- 
able limit for E" estimation accuracy. In this paper we 
show that using a simple line shape analysis — namely 
a Voigt fitting procedure — one can purify the intrinsic 
interaction effects from the extrinsic influence of the ex- 
perimental setup and, therefore, uncover true parameters 
of low energy part of quasiparticle spectrum. Applying 
the procedure to the nodal photoemission spectra from 
Bi 2 Sr 2 CaCu208+<5 (Bi-2212), we determine true values of 
the impurity scattering, as well as the energy and temper- 
ature dependence of the scattering rate in close vicinity 
to the Fermi level. 

We analyze the spectra from pure Bi-2212, lead-doped 
superstructure free Bi(Pb)-2212 M IE 0> as well as 
Bi-2212 doped with Zn and Ni Here we focus on 
the spectra measured along the nodal (it, tt) direction 
where the 5x1 superstructure is well resolved 01 an d at 
27 eV excitation energy at which the contribution from 
the bonding band is essentially suppressed -la- The ex- 
perimental details can be found elsewhere [SIS 0, IE E} • 

Fig.^, introduces the essentials of the nodal spectra 
analysis. The blurred region represents the "ARPES 
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FIG. 1: (a) Nodal spectra parameters: Bare band dispersion 
(dashed line) and renormalized dispersion (points) on top of 
the spectral weight of interacting electrons ("ARPES image"). 
Solid line represents a single MDC at ui, arrows indicate its 
FWHM and the self-energy parts derived for given MDC. (b) 
Two MDC's taken at -0.1 eV and -0.01 eV, and their fits to 
Lorentzian (solid lines) and Gaussian (dashed line) functions. 



image" — the photocurrent intensity, I(k,uj), which, over 
the occupied states, can be well approximated by the 
quasiparticle spectral function 



1 



A(k,Lu) = 



7r(w-£(/c)-£'(w)) 2 + £"H 2 ' 



(2) 



where e{k) is the bare band dispersion along the nodal 
direction. Since the ARPES image became a unit of in- 
formation in modern photocmission, the advantages of 
the analysis of the ARPES spectra in terms of the mo- 
mentum distribution curves, MDC = I(k,u> = const), 
had been immediately realized . The main advantage 
comes from the fact that A{k) has a simple Lorentzian 
lineshape as long as momentum dependence of the self- 
energy and bare Fermi velocity, vp — de{k)/dk, can be 
neglected 0. The latter, as well as Eq.||IJ, has been 
shown to be valid for the nodal direction of Bi-2212 up to 
0.3 eV binding energy by means of Kramers-Kronig self- 
consistency of the self-energy parts |5| . The relations of 
bare dispersion and self-energy parts with the parameters 
of an MDC at given u are shown in Fig.^,, though, in 
this paper we focus in E", which, at low binding energy, 
is simply proportional to the MDC width (half width at 
half maximum): E" = vpW. 

The analysis presented in Ref . |5| also allows to estimate 
contribution of the experimental resolution, although, 
due to a number of parameters involved, such an esti- 
mate is not very precise. Nevertheless, the close similar- 
ity between I(k,oj) and A(k, oj) at higher binding energies 
suggests the way to recover true quasiparticle spectrum 
also in close vicinity to Ep. It was noticed [jjj] that, 
when approaching the Fermi level, the lineshape of MDC 
measured along the nodal direction evolves from almost 
ideal Lorentzian to more Gaussian-like, which, evidently, 
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FIG. 2: Constituents of MDC width: (a) energy depen- 
dence of MDC width presented by half width at half max- 
imum (HWHM) of its fit to the Voight profile, W v , (Wi or is 
the result of simple Lorentzian fit), and its Lorentzian, Wl, 
and Gaussian, Wg, constituents; (b) variation of the above 
mentioned parameters (at ui = 0) with temperature. 



is a result of convolution of the photoemission signal with 
the total response function of the setup. Fig^illustrates 
such a Lorentzian to Gaussian crossover: While the MDC 
taken at lj = —0.1 eV is almost perfect Lorentzian, the 
MDC measured closer to the Fermi level at a; = —0.01 eV 
can be fitted to neither Lorentzian nor Gaussian but to 
a convolution of these two — the Voigt profile. 

Fig. |21 illustrates effectiveness of the fitting procedure 
for the nodal spectra analysis. The parameters of the fit 
are shown in Fig. E^l as a function of energy for Bi-2212 
OP89 sample measured above T c . Here Wv is the width 
of the Voigt profile, while Wl and Wq are the widths 
of its constituents, the Lorentzian and Gaussian, respec- 
tively. Wior represents the "traditional" MDC width used 
in previous data analysis — the width of the Lorentzian 
from the pure Lorentzian fit. Fitting MDC's to the Voigt 
profile instead of to the Lorentzian introduces one addi- 
tional parameter, Wg ■ In the fitting procedure we use the 
Voigt function 14] implemented in IGOR Pro (Wavemet- 
rics Inc.) which can be approximated by a simple relation 
between the above mentioned parameters: 



W T 

W V = V(W L ,W G ) = -± 



wl 



wl. 



(3) 
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FIG. 3: Energy dependence of the scattering rate, (a) Optimally doped Bi-2212 below T c : MDC width (HWHM) measured 
with 27 eV synchrotron radiation Q (filled squares) and the true width of the quasiparticle spectrum for the same sample 
purified from the resolution effect (filled triangles), to compare to the data for similar sample measured with 6 eV laser light 
\u\ (open squares), (b) Evolution of W(uj) with increasing temperature and overdoping. (c) The same data shown in a reduced 
dimensionless scale, dashed lines represent a parabolic fit to the low T data and a linear fit to the hight T data, respectively. 



Deviation of this approximation from real convolution is 
maximal when Wl ~ W G but is less than 1.2%. 

As it is expected for contribution of the experimen- 
tal resolution, the Gaussian width is essentially u>- 
independent in the actual range of interest, —0.1 eV < 
uj < eV: W G {uj) = W G . So, the energy dependence of 
Wy is accumulated in Wl, which represents now a true 
width of the quasiparticle spectral function. At higher 
binding energies (ui < —0.1 eV), Wg{u>) starts to increase 
that, as it follows from our simulations, can be explained 
by a non-linear bare dispersion. For even lower energies 
(< —0.2 eV), due to critical lowering of the signal to noise 
ratio, the fit becomes unstable in distinguishing the line- 
shape type. Therefore, in order to reduce the experimen- 
tal uncertainty we ht the Wg(uj) to Wq on [—0.1,0] eV 
energy range and define the true Lorentzian widths W[ 
from Eq. (1): W v = V(W[,W G ). In the following, we 
omit "/" but discuss exactly the W[(cu) function as the 
most careful representative of the true scattering rate. 

Fig. 13) represents not usual but interesting example 
that demonstrates the efficiency of the Voigt fitting pro- 
cedure. It shows the values WV(0), Wl(0), and Wg de- 
termined for the same sample at different temperatures. 
In this case, the APRES spectra have been measured 
during heating that implies a motion of the sample in 
respect to the beam spot due to thermal expansion of 
the manipulator. Evidently, this has resulted in different 
angular resolution in each point, that appeared as a ran- 
dom distribution of Wg(T) while the recovered Wl(T) 
function has found to be monotonic. 

We have applied the procedure to a number of ARPES 
spectra taken from Bi-2212 samples of different doping 
levels at different temperatures. In following we dis- 



cuss the energy dependence of the true scattering rate, 
Wl{uj), its zero energy value, Wl{0), and its dependence 
on doping and temperature. 

In respect to Wl(uj) problem, Eq. helps to make 
an important remark. In a low scattering limit, when 
W L < W G , W v ~ W G + W L /2. In other words, if, for 
example, Wl{<+>) oc uj, such a linear dependence cannot 
be camouflaged by the resolution. So, the use of the 
Voigt profile instead of the Lorentzian does not influence 
the qualitative result for the dependence of the scatter- 
ing rate on binding energy — one can see no principal 
difference between the Wj,(w) and W\ a ^{uj) functions in 
Fig- [2^ — but becomes highly important if we are inter- 
ested in its careful quantitative analysis []|. 

With Fig. |3 we discuss the energy dependence of the 
scattering rate. In panel (a) we compare the MDC width 
for OP89 sample, measured below T c with 27 eV syn- 
chrotron radiation , to the data for the near-optimally 
doped Bi-2212, measured with 6 eV laser light |Tjj . Be- 
sides different offsets, two datasets look very similar, [see 
also panel (c)] . Moreover, when the width of the quasi- 
particle spectrum for OP89 sample is purified from the 
resolution effect, it almost coincides with the laser data. 
This supports the conclusions about validity of the sud- 
den approximation and that the ARPES spectra on Bi- 
2212 represent its bulk properties. 

To reveal the asymptotic behavior of the scattering 
rate at low energy, in panel (c) we replot these data in 
form of a"((j) = [£"(w) - S"(0)]/cj. To get E" we just 
multiply Wl by vf = 4 eVA Q. One can see that, below 
T c at low energy, <j"{oj) has a non- vanishing offset and 
a quadratic term. In S"(w), these terms correspond to 
~ uj and ~ uj 3 terms, respectively. At higher temper- 



4 



OD75 




50 100 150 200 250 300 
T (K) 

FIG. 4: Temperature dependence of the scattering rate. Er- 
ror bars are determined by mean square deviation of Wq, 
solid lines show the results of 2nd order polynomial fitting. 

atures and hole doping [see panels (b) and (c)], E"(u;) 
becomes purely quadratic. At very low energy (|o;| < 
50 meV), both terms below T c can be explained by in- 
fluence of the d-wave superconducting gap which causes 
the density of states (DOS) to become linear at low en- 
ergy: the linear term is expected for the elastic forward 
scattering and the cubic term can be associated with in- 
elastic spin-fluctuation scattering Also, in panel (c), 
we draw the "quasiparticle limiting line" to show how, in 
fact, well defined are the nodal quasiparticles in Bi-2212. 

Fig. 01 shows Wl(0,T) dependences for several sam- 
ples. It is important to stress that the scattering rate, 
even when purified from resolution, does not vanish at 
lj -> and T -> 0. The residual, Eg = v F W L (0, 0), can 
be naturally explaine d by scattering on impurities. It has 
been shown recently that such a high value (Eg > 
16meVor W L (0,0) > 4 X 1CT 3 A" 1 ) can be caused by a 
forward scattering on out-of-plane impurities rather than 
by usual unitary scattering on in-plane impurities. How- 
ever, we should note, that the determined Eg values are 
more than 7 times lower than the values compared ear- 
lier to resistivity measurements , therefore the actual 
discrepancy between transport and ARPES experiments 
is much lower. 

Another conclusion that can be derived from the off- 
set values in Fig. 0] is that Eg does not show a system- 
atic dependence on doping level (OD75, OP89 and UD78 
samples), but is rather sensitive to impurity concentra- 
tion (here UD76 represents a sample, highly underdoped 
by annealing in He, which removes oxygen from SrO 
planes introducing a number of out-of-plan defects, while 
OP86 sample doped with 1% of Zn gives an example 



of the in-plane impurity). This allows to conclude that 
the increase of the residual resistivity with underdoping 
still cannot be explained by an increase of the scattering 
rate but only by a decrease of the charge carrier density, 
which is inconsistent with "large" Fermi surface with area 
~ 1 — x , providing therefore an indirect evidence for 
the phase separation in underdoped cuprates. 

T-dependence of E" deserves a separate investigation 
and discussion. One of the surprising results is that, 
within the accuracy ~ 10 meV, there is no drop of E" 
detected at T c . This can be a consequence of the pseudo- 
gap, which smoothes the T-dependence of the phase 
space available for scattering. Such T-dependence is a 
reason for aw-T asymmetry of E"(w, T), and compari- 
son of E"(0, T) and T,"(ui, 0) can be used to determine the 
phase space evolution. It is also important to understand 
the reason for different dependences presented in Fig. 01 
While for the most samples the variation in E"(0,T) 
can be explained by different scattering on impurities, 
the much stronger T-dependence for the Y-doped sample 
needs further investigations. Answering these questions 
require data of much better accuracy. The described pro- 
cedure eliminates the sample-related uncertainty, so, the 
final accuracy depends only on experimental statistics. 
Therefore, only time is needed to clarify these points. 
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